Effects of momentum conservation on the analysis of anisotropic flow 
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We present a general method for taking into account correlations due to momentum conservation 
in the analysis of anisotropic flow, either by using the two-particle correlation method or the standard 
flow vector method. In the latter, the correlation between the particle and the flow vector is 
either corrected through a redefinition (shift) of the flow vector, or subtracted explicitly from the 
observed flow coefficient. In addition, momentum conservation contributes to the reaction plane 
resolution. Momentum conservation mostly affects the first harmonic in azimuthal distributions, 
i.e., directed flow. It also modifies higher harmonics, for instance elliptic flow, when they are 
measured with respect to a first harmonic event plane such as one determined with the standard 
transverse momentum method. Our method is illustrated by application to NA49 data on pion 
directed flow. 

PACS numbers: 25.75.Ld 25.75.Gz 



I. INTRODUCTION 

In heavy ion collisions, the transverse momenta of out- 
going particles are correlated with the direction of the 
impact parameter (reaction plane) of the two incoming 
nuclei. For example, at energies above 100 MeV per nu- 
cleon, nucleons in the projectile rapidity region are de- 
flected away from the target. This collective effect is 
usually referred to as directed (or sidewards) flow. It is 
characterized by the first Fourier coefficient of the parti- 
cle azimuthal distribution, v\ |Q : 



(1) 



where <f> denotes the azimuthal angle of an outgoing par- 
ticle, &rp is the azimuth of the reaction plane, and an- 
gular brackets denote an average over events and over all 
particles in a given transverse momentum and rapidity 
window. Differential flow is characterized by v\(pT,y) 
that describes the directed flow in a narrow transverse 
momentum (pt) and rapidity (y) interval. 

The most common observable to quantify directed 
flow, introduced by Danielewicz and Odyniec in 1985 
is the mean transverse momentum projected on the re- 
action plane, (p x ) , as a function of rapidity. In terms 
of vi, it can be written as (p x ) = (pt vi(j>T,y)), where 
the average in the right-hand side is taken over px- This 
observable is believed to be sensitive to the nuclear com- 
pressibility and to in- medium cross sections |3j. Experi- 
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mentally, (p x ) is extracted from azimuthal correlations 
between the outgoing particles, under the assumption 
that these correlations are only due to flow. However, 
it was soon realized that global momentum conservation 
also induces azimuthal correlations. Different ways to 
correct for this effect were devised in Refs. ||, ||. These 
methods are still in use at intermediate |(| [7J as well as 
at relativistic energies [^], [| Generally, corrections 
are important for weak flow and/or small systems, when 
there are physical or detector differences between the for- 
ward and backward hemispheres in a reaction. 

More detailed information can be obtained on the 
physics involved in the collision by studying the trans- 
verse momentum dependence of directed flow |ll[, v\(pt) 
in a given rapidity window. Detailed theoretical stud- 
ies have been presented for this differential directed flow 
]l2| , |l3| , fl4| . This was accompanied by new analyses of 
experimental data, and a wealth of results are now avail- 
able from FOPI at SIS E877 at AGS fTfJ and NA49 
at CERN p7y . In these analyses, unfortunately, correla- 
tions arising from momentum conservation are not taken 
into account. Recently, it was shown that they may 
in fact be of the same order of magnitude as correlations 
due to flow, leading to a reevaluation of the "standard" 
flow analysis results. 

In this paper, we show how to modify standard flow 
analysis techniques in order to take into account momen- 
tum conservation. Two flow methods are used experi- 
mentally. The first one relies on a study of two-particle 
correlations Q. This simple case is discussed in Sec. Q. 
The second method, which is by far the most common, 
involves the construction of a "flow vector" used as an es- 
timate of the reaction plane ^ . In order to measure 
the flow of a particle, one correlates it to the estimated 
event plane of the flow vector, and corrects this correla- 
tion by dividing by a "resolution" factor which takes into 
account the uncertainty of the reaction plane determina- 
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tion. This procedure is recalled in Sec. LEI. In Sec. IV, 
we show how this method can be modified to take into 
account momentum conservation. 

The momentum conservation leads to two effects, 
whose magnitude is controlled by a parameter / which is 
roughly the square root of the fraction of all particles used 
in the estimate of the event plane. First, momentum con- 
servation contributes to the correlation between a parti- 
cle a nd the flow vector, proportionally to the parameter / 
(Sec. IV B). This spurious correlation can be elimin ated 
through a redefinition of the flow vector (Sec. IV A ). In 
addition, momentum conservation affects the resolution 
of the event plane, a lthou gh the effect is smaller, as it is 
quadratic in / (Sec. IV C| ). This second effect may also 
bias higher harmonics measured with respect to the first 
harmonic event plane. (Note, however, that momentum 
conservation has no effect when elliptic flow is measured 
with respect to the second harmonic event plane.) Our 
results are discussed in Sec. |v|; in particular, we show 
that the correlation from momentum conservation van- 
ishes if the detector acceptance is symmetric with respect 
to midrapidity, because particles belonging to the for- 
ward and backward he misp heres are given weights with 
opposite signs. In Sec. [Vi], the method is illustrated by 
NA49 data on pion directed flow. Finally, technical de- 
tails are exposed in the Appendixes |X| and [b]. 



II. TWO-PARTICLE CORRELATION 
TECHNIQUE 

Since the orientation of the reaction plane is not 
known, v\ is extracted from azimuthal correlations be- 
tween outgoing particles. The simplest way consists in 
using two-particle azimuthal correlations p9| . These are 
related to v\ by the identity: 

(cos(0i - (j> 2 )) = (cos(0i - $ RP ) cos(0 2 - §rp)) 

= (cos(0l - <f>R P )) (C0s((f> 2 - <f>R P )) 

= »i(l)«!(2), (2) 

where uj,(l) and t>i(2) denote the flow coefficients asso- 
ciated with each particle. For spherical nuclei, symme- 
try with respect to the reaction plane ensures that the 
(sin(0 — $rp)) terms are zero. In going from the first to 
the second line, we have assumed that all azimuthal cor- 
relations are due to flow or, equivalently, that azimuthal 
angles with respect to the reaction plane <fii — $pp and 
4>2 — &rp are statistically independent. Recent imple- 
mentations of the two-particle correlation method can 
be found in Refs. § ||. 

As recalled in the introduction, however, there are also 
azimuthal correlations induced by global momentum con- 
servation, so that the measured correlation is in fact the 
sum of two terms(when both effects are small) 



where the last term, due to momentum conservation, is 
evaluated in Appendix A of Ref . [|18| : 



(cos 



(01 - <h)) 



mom. cons. 



ptipt 2 
~N(&)" 
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where N denotes the total number of particles produced 
in the collision, and the average value (p?) 1S similarly 
taken over all produced particles. Since one usually de- 
tects only a fraction of the produced particles, the denom- 
inator of Eq. ([|) should be estimated using a model, or 
by an extrapolation of available data to full phase space. 

In order to take into account the correlations due to 
momentum conservation in the analysis of directed flow, 
one simply subtracts the contribution due to momentum 
conservation, Eq. (^), from the measured two-particle 
correlation Eq. (||), so as to isolate the correlation due 
to flow. This is the procedure followed in Ref. Jll| to 
estimate the effects of momentum conservation on the 
analysis of directed flow at SPS. 



III. FLOW- VECTOR METHOD 

Most analyses of directed flow rely on the construction 
of the event flow- vector 10, |2Cf| defined as: 



JV 

Q = >'/(', U 



j 3 ' 



(5) 



(cos(0i - 02 )) = t>i(l)t>i(2) + (cos(0a - 02 )) 



mom. cons. 1 

(3) 



where the sum runs over all produced particles, Uj de- 
notes the unit vector of the particle transverse momen- 
tum, Uj = e x cos 4>j + e y sin (fij = p^j /pTj , and Wj is 
a weight depending on the particle type, its transverse 
momentum ppj and rapidity y,j. The weights, Wj, are 
usually given non zero values only in regions of good de- 
tector acceptance. 

The azimuthal angle of the flow vector, $, is consid- 
ered an estimate of the reaction plane azimuth, $>pp. 
Therefore, the azimuthal angle of Q with respect to the 
reaction plane is usually referred to as the uncertainty in 
the reaction plane determination, and will be denoted by 
A$ = $ — $7?p. One then writes an equation similar to 
Eq. (J^) for the azimuthal correlation between a particle 
with azimuth <p and the flow vector: 

(cos(0 - $)) = (cos(0 - <$> RP - $ + <& RP )) 
= (cos(0 - $_rp)) (cos A$) 
= vi (cosA$) . (6) 

This equation once again relies on the assumption that 
azimuthal correlations are only due to flow. To avoid 
autocorrelation one should subtract the contribution of 
the particle under study from the flow vector. Since the 
flow vector, Eq. (|^), involves a summation over many 
particles, the correlation between a particle and the flow 
vector, Eq. (^|), is usually much stronger than the corre- 
lation between two particles, Eq. (0). This is one of the 
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main reasons why the flow vector method is more com- 
monly used. However, correlations due to momentum 
conservation are also stronger, so that corrections due to 
momentum conservation are of the same magnitude with 
both methods, as we shall see in Sec. In the rest of 
this paper <&rp will be taken as zero; that is, the reaction 
plane is along the a;- axis. 

The next step in the analysis is to evaluate the "reso- 
lution" (cos A$). This factor results from a competition 
between flow, which tends to align Q in the direction 
of the reaction plane, and statistical fluctuations, whose 
relative magnitude decreases like 1/y/N. A quantitative 
estimate can be obtained with help of the central limit 
theorem. This theorem allows one to write the distribu- 
tion of the flow vector as a Gaussian distribution in the 
coordinate system where the orientation of the reaction 
plane is fixed |TJ p2[ : 



dP 
dQ 



7T(T 



W (Q ~ <Q>)2 



a- 



(7) 



The average value (Q) is aligned with the reaction plane, 
i.e., with the £-axis, while the dispersion a is due to 
statistical fluctuations. Writing Q = Q(e x cosA$ + 
e v sin A<fr), and integrating over Q and A$, one ob- 
tains lEoLpll: 



( C0S A$) = ^ X e-* 2/2 [/o(x 2 /2) +/i(x 2 /2)] , 



(8) 



where Iq, I\ are modified Bessel functions and x is the 
dimensionless parameter x = (Q) ' e x/c, which char- 
acterizes the relative magnitude of flow and statistical 
fluctuations. In what follows, x wm b e referred to as the 
"resolution parameter" (34| . 

This parameter x is usually estimated from the corre- 
lation between two "subevents" || : the set of N particles 
is divided randomly into two equivalent subsets of N/2 
particles, labeled a and b. One then constructs two flow 
vectors Q a and Qb defined as in Eq. (||) , with the summa- 
tion running over the corresponding subset. Flow induces 
an azimuthal correlation between Q a and Qb, since both 
are correlated with the reaction plane. Then, when all 
azimuthal correlations are due to flow, the distribution of 
$ a — depends only on the resolution parameter x- To 
find the functional form of this relationship, one writes 
the probability distribution of Q a in a form similar to 
Eq. (0): 



dP 



i 



■ exp 



(Qa - (Qa))' 



(9) 



and a similar formula for Qb with (Qb) = (Qa) = (Qsub) 
and <ib = & a = Osub since the subevents are equivalent. 
The distribution of the relative angle dP/d($ a — $£,) de- 
pends only on the resolution parameter associated with 
either subevent Xsub = (Qsub) ■ e x /cr su b. Its analytic 
expression can be found in Refs . p3| , pi[ , and provides a 
perfect fit to experimental data ]25|. The parameter Xsub 



is usually obtained from the square root of the average 
cosine of the relative angle, y/ (cos($ a — which is 

simply the resolution of a subevent plane, and is given 
by Eq. (||) with x replaced by Xsub- Equation (H) must 
then be inverted in order to obtain Xsub [ p0| ■ In a simpler 
way, Xsub can be obtained directly from the fraction of 
events with the relative angle between subevents |<& a — $b| 
larger than tt/2 [|| M: 



P(|*„-$ 6 | >7T/2) 



1 



"Xsub 



(10) 



This method is used for example in Ref. (25|. If Xsub 
is significantly larger than unity, however, only a small 
fraction of events have |$ — <E>;,| larger than tt/2, and 
the resulting statistical error on Xsub is larger than with 
the previous method. In the limiting case Xsub 3> 1, 
the distribution of the relative angle $ Q — $6 is a Gaus- 
sian, and Xsub can be directly obtained from its width, 

<(<fa-$ fc ) 2 ) = l/Xs 2 ub- 

The last step is to relate the parameter Xsub, corre- 
sponding to one subevent, to x, corresponding to the 
whole event. Since each subevent contains half of the 
particles, (Q su b) = (Q) /2. If correlations between parti- 
cles are only due to flow, a similar scaling occurs in the 
statistical fluctuations: of ub = er 2 /2, so that 



X = Xt 



ubv2- 



(11) 




The flow- vector method is defined by Eqs. (||), (| 
and a technique for obtaining Xsub- Equations (f 
and (|ll|), must be modified when taking into account 
momentum conservation. This is studied in Sec. IV. 

The first harmonic flow vector, Eq. (||), can also be 
used to measure higher harmonics of the azimuthal dis- 
tribution v n — (cos n<j)) for n > 1, with the help of an 
equation similar to Eq. (||): 



(cosn 



,_$)) 



(cosnA$) 



(12) 



The corresponding resolution (cos nA$) depends on the 
same parameter x |2(], p3l. T his method has been used 
to analyze elliptic flow V2 ]2q| relative to a first harmonic 
event plane. Since, as we have mentioned above, momen- 
tum conservation modifies the value of x, it must also be 
taken into account in this case. We come back to this 
discussion in Sec. 0. 



IV. MODIFICATIONS DUE TO MOMENTUM 
CONSERVATION 

Global momentum conservation induces correlations 
between outgoing particles and the flow vector. In the 
same way as Eq. (||) is replaced by Eq. (||), Eq. (0) be- 



(cos(0 - $)) = vi (cos A$) + (cos((/> - $)) 



mom. cons. 



(13) 



4 



There are two methods to take into account the additive 
term. The first one consists in redefining the Q-vector 
[i.e., modify $ in Eq. Jl3|)] in each event, in such a wa y 
that the additional term no longer appears (Sec. IV A ). 
Alternatively, one can compute explicitly the additive 
term and subtract it in order to isolate the correlation 



due to flow (Sec. IV B). Furthermore, momentum conser- 
vation affects the resolution (cosA$), and the methods 
used to est imate the resolution must be modified accord- 
ingly (Sec. HVCj ). 



A. Shifted flow-vector method 

The correlation due to momentum conservation be- 
tween the flow vector and a particle with momentum pt 
can be characterized by the average value of u • Q, where 
u = Pt/pt- Using the definition of the Q-vector, Eq. (||), 
and the expression of the two-particle correlation due to 
momentum conservation, Eq. (|4|), one obtains: 



(u-Q) 



mom. cons. 



-Pt ■ 



[WPT) 

(Pt) 



(14) 



where averages in the right-hand side run over all the 
produced particles. This suggests one could redefine the 
flow vector as 



<4 — <4 + , 2 r Pt, 

(Pt) 



(15) 



so that momentum conservation does not contribute to 
the correlation between Q' and the particle under study. 
As a matter of fact, when there is no flow the shifted flow 
vector Q' and px are now uncorrelated by momentum 
conservation: any correlation between Q' and px arises 
from flow, because here and throughout the paper, we 
neglect other sources of nonflow correlations. 

A similar method was proposed in Ref. Q , in the con- 
text of intermediate-energy collisions, where the shifted 
flow vector was defined as 



14 — 14 H Pt, 



(16) 



where mj is the mass of particle j. If (p^) of different 
particles scales with their mass m, our definition Eq. Jl5| ) 
reduces approximately to Eq. (|l6|). That is the case in 
particular for a thermalized nonrelativistic gas at rest 
with (almost) no flow, as in Ref. §. However, Eq. @ 
is no longer valid when there is flow, and/or relativistic 
effects, whereas Eq. (fT5J) still holds. 



B. Explicit subtraction 

Instead of modifying the flow vector event-by-event, 
one can as well keep the original definition, Eq. (^), and 
subtract the contribution of momentum conservation to 



the measured correlations, Eq. (fl3|), at the end of the 
analysis. 

Let us first estimate the order of magnitude of the last 
term of Eq. (|l3|). We assume for simplicity that the par- 
ticle under study has no flow, v\ = 0, so that only the 
correlation due to momentum conservation remains. We 
then make the following approximation: 



(cos(</> - <&)) 



mom. cons. 



Q 



Q 



(u-Q) 

(Q) 



(17) 



with u = px/pt- The average dot product is given by 
Eq. (|l4[) . Furthermore, the average length of the flow 
vector (Q) is approximately equal to the r.m.s. length 
yj (Q 2 )- Neglecting correlations between momenta, one 
obtains from the definition, Eq. (|), (Q 2 ) ~ N (w 2 ). All 
in all, the order of magnitude of the correlation due to 
momentum conservation, Eq. (|17|), reads 



d>)) 



Pt 



mom. cons. 



where / is a dimensionless quantity given by 
(wp T ) 



/, (18) 



/ = 



y/WTW) 

(wPt) c 



M 



{w 2 ) Q N{p 2 T Y 



(19) 



and the subscript Q refers to those M particles used for 
the Q-vector (for which weight Wi 7^ 0). The weight w 
vanishes for particles not seen in the detector, so that 
(w 2 ) and (wpr) scale like the fraction of detected parti- 
cles. Hence, / scales like the square root of the fraction 
of detected particles. This is seen more clearly in the 
second line of Eq. (|l9|) where there is a factor yjM/N. 
If to is taken as +1 for particles detected in the for- 
ward hemisphere and —1 for the backward hemisphere, 
then (w PT ) Q = (M F {p T ) F - M B (p T ) B )/M, where F 
and B refer to particles used for the Q-vector from the 
forward and backward hemispheres, respectively, so that 
M = Mp + Mb- In this case the quantity (w 2 ) n = 1. 

The estimate in Eq. ( |l8| ) is only an order of magnitude. 
A more accurate expression can be obtained under very 
general assumptions. We show in Appendix [A| that the 
actual value of the correlation due to momentum conser- 
vation between a particle and the flow vector is 



$)) 



mom. cons. 



Pt 



f 



: e-* 2 /W/2), (20) 



where \ = (Q) ' e x/c has been defined in Sec. |HI|. It 
is worth noting that the Eq. ( p0| ) result now depends 
on x, expressing the effect of flow in reducing the cor- 
relation due to momentum conservation. However, the 
dependence is weak for small x (i- e -, m the case of poor 
resolution) . 



5 



C. Correction to the reaction plane resolution 

We now discuss the influence of momentum conserva- 
tion on the determination of the resolution parameter \. 
The corrections discussed below hold for the method of 



Sec. IV A as well as for the method of Sec. 1VB 



As recalled in Sec. Ill, the resolution is determined 
from the correlation between subevents, i.e., from the 
distribution of the relative angle $ a — between the 
subevents. We show in Appendix [b| that momentum con- 
servation contributes to this correlation such that the dis- 
tribution of <& a — <&b now depends both on the resolution 
parameter of the subevent Xsub, and on the parameter / 
introduced in Eq. (|l9|), in a non trivial way: taking into 
account momentum conservation does not only amount 
to changing the value of Xsub- This is because the shape 
of the distribution of the relative angle <i> a — $t depends 
on whether the correlation is due to flow or to other ef- 
fects. This is illustrated in Ref. ^7j, where the distribu- 
tion of $ a — $h is plotted in the two limiting cases when 
the correlations are only due to flow (/ = 0), or only to 
other, nonflow effects (xsuh — 0). 

As recalled in Sec. 
extracting Xs 



III 



there are several methods of 
ah from the distribution of $ a — The 
usual method relies on the calculation of (cos( ( I> a — $&)). 
If / is much smaller than unity, Eq. (||) becomes: 

(cos($ a -$ fc )) 1/2 = 



Xsub l/ofrsub/ 2 ) 

^1 (Xsub 



y [^(Xsub/2) 



h (Xsub 

1/2 

/2)] 



2 /2)] 2 



(21) 



Once the left-hand side and / are known, Xsub can be 
obtained by solving this equation. This expression holds 
only to leading order in /. We were not able to obtain 
an analytic formula for arbitrary /. 

It may be simpler to extract x su b from the fraction of 
events for which |$ a — <J>(,| > ir/2. Then, the correction 
for momentum conservation can be implemented using 
the following exact formula, derived in Appendix [b|: 

P(|$ Q _$ b |>7r/2) = ^^e-*»>>, (22) 

instead of Eq. (|To|). If the resolution parameter is large, 
Xsub > 1) this method leads to larger statistical errors 

If Xsub 



1, the distribution of 



as recalled in Sec. III. 
<I> Q — is Gaussian and Xsub can be obtained from its 
width: 

((3>a-<& fc ) 2 )= 2 2 ■ (23) 

A su bV z J I 

Finally, the relation between the resolution parameter 
of the subevent, Xsub: and that of the whole event, x, 
is easily obtained using Eqs. (|A9|), ([B2|), and the scaling 
relation (Q sub ) = (Q) /2: 



X — Xsubl 



' 2-/ 2 
I-/ 2 



(24) 



When correlations are weak (i.e., when both / and x 
are small compared to unity), the event plane resolution 
can be corrected for momentum conservation through the 
following simple formula: 



X 2 (corrected) = x 2 (uncorrected) + / 



(25) 



This approximation can be used instead of Eqs. ( 21 )/ fl22] ) 
and {vA ) on the x of the full event plane as obtained in 



Sec. Ill 



Let us summarize the results obtained in this Section. 
The crucial parameter, which determines the magnitude 
of corrections due to momentum conservation, is the di- 
mensionless parameter /, defined in Eq. (|l9|). Once this 
quantity has been estimated, the subevent plane res- 
olution parameter, Xsub, can be obtained from either 
(cos($ a — $;,)) using Eq. if / <C 1, or from the frac- 
tion of subevent angle differences greater than it/ 2 using 
Eq. (p2|). Or, if the resolution parameter is very large, 
from Eq. (^3|). The full event plane resolution parame- 
ter, x, is obtained from Eq. (p4|), and the full event plane 
resolution from Eq. (|^). Finally, subtracting from the 
measured values of (cos(0 — <&)) the contribution due to 
momentum conservation, Eq. (|20|), and dividing the re- 
sult by the event plane resolution, one obtains the final 
i>i values. When / vanishes, one recovers the formulas of 
the standard flow analysis, recalled in Sec. III. 



V. DISCUSSION 

Corrections arising from momentum conservation are 
twofold. First, they contribute to the correlation between 
a particle and the event plane. This contribution, which 
increases linearly with the particle transverse momen- 
tum must be subtracted from the measured corre- 
lation so as to isolate the correlation due to flow; see 
Eq. (jl^). Next, momentum conservation affects the reso- 
lution of the event plane, which results in a global multi- 
plicative factor. Both effects must be taken into account 
in the analysis of directed flow, but only the second one in 
the case of higher harmonics relative to a first harmonic 
event plane. 

All effects of momentum conservation disappear if the 
parameter / defined in Eq. ( |l9| ) vanishes, which occurs 
if (wpx) = 0. Since w is usually an odd function of the 
center-of-mass rapidity (wpx) vanishes as soon as the 
detector acceptance is symmetric with respect to midra- 
pidity, as was already mentioned in Refs. ^0|. Another 
theoretically interesting limit is the case when all parti- 
cles are detected and, in particular, when w = px, which 
implies / = 1. In this case, the flow vector, Eq. (||), is 
the sum of all transverse momenta, which vanishes due 
to momentum conservation. This is why the corrections 
derived in the previous Section diverge when / goes to 
unity. 

When the weight keeps a constant sign, the correction 
due to momentum conservation is of the same relative 
magnitude within the flow vector method [Eq. (|l^)] as 



6 



within the two-particle method [Eq. (0)]. In the latter 
method, the correlation due to flow is of order v\ while 
the correlation due to momentum conservation is of or- 
der 1/N. The flow vector method enhances both terms 
by a factor of order y/M, where M is the number of par- 
ticles within the detector acceptance, which are used in 
constructing the flow vector Q. Roughly speaking, this 
is because the flow vector is analogous to a random walk 
of M steps: the length of the walk is approximately \[M 
times the size of a step, and similarly the correlation be- 
tween Q and a given particle is \f~M times the correlation 
between two particles. This can be seen more explicitly 
from the formulas derived above. We first estimate the 
correlation due to flow, i.e., the first term in the right- 
hand side of Eq. ( fi"3| ) . Since (cos A<£>) is proportional to \ 
for small \ ( see Eq. ([§) and the figures in Refs. |^0|, p3|), 
and x ~ v i VM, this term is indeed of order v\ \/M. Let 
us now estimate the second term, whose order of mag- 
nitude is given by Eq. ([l8]). We choose for simplicity a 
weight w — pt within the detector acceptance, and zero 
outside. The parameter / in Eq. ( |l9| ) then reduces to 
/ ~ y/ M/N. Introducing this value in the estimate Eq. 
(|l8|), we find a correlation of order \/~M/N, as expected. 

The correction to the reaction plane resolution, on the 
other hand, is generally of order f 2 , which is roughly the 
fraction of the total number of particles used in estimat- 
ing the reaction plane. Note that the correction to the 
resolution is not the same for all harmonics of the az- 
imuthal distribution. Indeed, the resolution (cosnA$) is 
more sensitive to \ when n > 1, so that effects are larger 
for higher harmonics relative to the first harmonic event 
plane. 

Let us now discuss how the corrections vary with the 
centrality of the collision. The contribution of momen- 
tum conservation to the correlation between a particle 
and the flow vector scales with the number of particles N 
like 1 /N. Therefore, momentum conservation may in fact 
dominate the measured correlation for the most periph- 
eral collisions. For central collisions, the flow vanishes 
by symmetry and only momentum conservation remains, 
but it is likely to be smaller due to the larger value of 
N. In contrast, the parameter /, which determines the 
correction to the event plane resolution, is to a good ap- 
proximation independent of centrality. This is because 
the fraction of particles in the detector acceptance varies 
little with centrality. One could use this scaling prop- 
erty to estimate the contribution of momentum conser- 
vation (and other nonflow effects) experimentally. This 
property was used in Ref. [p8| , where the contribution of 
momentum conservation was estimated from the correla- 
tions between four subevents. 

Besides analysis of the centrality dependence, one can 
observe the contribution due to momentum conservation 
by checking if v\ (pt ,y) = at midrapidity for all trans- 
verse momenta. In principle, the correlation of particles 
at midrapidity with other particles can be used to esti- 
mate the momentum conservation effect. Checking the 
dependence on pt shown in Eq. (H) , could also be useful. 



It is important to notice that there still exist other non- 
flow effects , as for instance correlations due to reso- 
nance decays or quantum effects. One can try to avoid 
them by performing appropriate cuts in phase space or 
using carefully chosen subevents j2^|. Alternatively, one 
can take them into account by performing a similar cor- 
rection to what we propose here, through a detailed mod- 
eling of nonflow effects jl8|, |3(j). Finally, if statistics 
are large enough, one can construct correlations between 
more than two particles: then, nonflow effects can be 
eliminated in a model independent way . 

VI. APPLICATION TO NA49 DATA 

To experimentally implement the correction we derived 
in the previous sections, one should estimate the various 
coefficients which enter Eq. (^0|). Since all emitted par- 
ticles are not detected, one needs a parameterization of 
the cross section for each type of emitted particle, and 
for each centrality bin, in order to compute the average 
(Pt) for all particles and the total multiplicity N. One 
can then check that the results make sense by plotting 
the rapidity dependence of the corrected directed flow, 
v\{y): if the signal which remains after correction is in- 
deed flow, it should cross zero at midrapidity, since v% is 
expected to be an odd function of y — ycu- 

FIG. 1: NA49 results on directed flow as a function of ra- 
pidity for charged pions from minimum bias 158^1 GeV Pb + 
Pb. Shown are vi before (squares) and after (circles) correc- 
tion for momentum conservation. Solid lines are polynomial 
fits. The open points have been reflected about midrapidity. 



4 



: \Nv 

^Ij^Tjri 


no corr. -. 
^V^. : 


: p, < 2 GeV/c 

: . ... i .... i .... i ... 


\*\ : 

with corr/Kj-J^ 
i .... i .... i .... i ... . : 



1 1.5 2 2.5 3 3.5 4 4.5 5 

rapidity 

As an illustration, the method was applied to the NA49 
data in Pb+Pb collisions at 158 GeV per nucleon [j32| . 
In the analysis, the event plane was determined from 
charged pions with px < 1-0 GeV/c, except for cen- 
trality bins 1 and 2, where the upper cuts were at 0.3 
and 0.6 GeV/c, respectively. The y values were be- 
tween 4.0 and 6.0. The particles were assigned weights 
w = \y — ycu\, always positive since all used particles be- 
long to the forward hemisphere. Then pions with known 
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Pt and y values were correlated with this event plane (be- 
ing careful to remove autocorrelations) to provide a ma- 
trix of observed v\ (px ,y)- The Xsub values were obtained 
from the correlation of subevent planes using Eqs. (|is| ) 
and (|^). In each (px, y) bin the momentum conservation 
correction was subtracted using Eq. (po|). The full event 
plane resolution was obtained from Eq. (g) and, accord- 
ing to Eq. (jl3]), the whole matrix was corrected for the 
event plane resolution to produce the double differential 
V\ values. The procedure was done for each centrality 
bin separately and then the minimum bias results were 
produced by taking the weighted average over centrality, 
weighting with the cross section. The cross sections as a 
function of pt, y, and centrality, had already been eval- 
uated and parameterized [j33| . Figure was produced by 
projecting this final matrix onto the rapidity axis, by tak- 
ing the weighted average over px, again weighting with 
the cross section. 



TABLE I: Listed for the six centralities are: N, the total mul- 
tiplicity over all phase space; (pr), the mean pi over all phase 
space in units of GeV 2 /c 2 ; M, the multiplicity of particles 
used for the event plane; /, the fraction defined by Eq. (|l9|); 
X, the resolution parameter for the full event plane; the res- 
olution of the full event plane; and the percent increase in 
the resolution due to momentum conservation. The first two 
centralities had more restrictive pr cuts. 



Centrality 


1 


2 


3 


4 


5 


6 


N 


2402 


1971 


1471 


1028 


717 


457 


(p 2 T ), GeV 2 / C 2 


0.32 


0.32 


0.31 


0.30 


0.28 


0.27 


M 


135 


184 


152 


108 


75 


46 


f 


0.07 


0.14 


0.17 


0.17 


0.17 


0.17 


X 


0.20 


0.25 


0.39 


0.45 


0.48 


0.47 


resolution 


0.18 


0.22 


0.33 


0.38 


0.40 


0.39 


res. increase, % 


6.3 


19.7 


10.9 


8.0 


7.1 


7.3 



The parameterized cross sections were also used to de- 
rive the (pt) f° r & H particles in Eqs. (19) and (po|). More 
specifically, in each centrality bin the average was ob- 
tained by taking 3/2 the sum of p\ of the charged pions, 
plus twice the sum for the charged kaons, plus the frac- 
tion of nucleons which are thought to be participants, and 
then dividing by the total multiplicity. The values of the 
parameters used in this momentum conservation correc- 
tion are shown in Table |. The (J^Px) = N (pr) values 
range from 773 to 122 GeV 2 /c 2 ; a slightly higher value 
was used in Ref. p8| , resulting in a somewhat smaller 
correction. The last line of the table shows that the in- 
crease in resolution due to momentum conservation is 
always less than 20%. 



However, Fig. [j] shows that the correction for the effect 
of momentum conservation on v\ is indeed significant, 
roughly an absolute value of 1% in the whole rapidity 
range under study. The corrected correlation does give a 
directed flow vanishing at midrapidity within error bars, 
as it should. In fact, if one did not know N and (p T ), one 
could force the data to go through zero at midrapidity as 
a way to evaluate this correction. 



Acknowledgments 

We would like to thank the NA49 Collaboration for 
permission to use their data and, in particular, Alexan- 
der Wetzler for performing the data analysis and Glen 
Cooper for parameterizing his cross section values. We 
also thank the anonymous referee for suggesting us the 
comparison with the method of Ref. N. B. ac- 

knowledges support from the "Actions deRecherche Con- 
certees" of "Communaute Francaise de Belgique" and 
IISN-Belgium. 



APPENDIX A: CORRELATION BETWEEN ONE 
PARTICLE AND THE EVENT FLOW VECTOR 

In this Appendix, we compute the azimuthal correla- 
tion between a particle with momentum px and the flow 
vector Q when both flow and momentum conservation 
are taken into account. Flow, which is defined as a cor- 
relation of outgoing momenta with the reaction plane, 
means here a correlation with the x-axis: the transverse 
momentum vector px in a given y and px window has 
a non-zero average value (pt) which is parallel to the 
x-axis, and the same holds for the flow vector Q. If the 
only correlation was due to flow, we could write the nor- 
malized iV-particle momentum distribution as a product 
of single-particle distributions: 



dP 



dpTi ■ ■ -dp TN 



Fi(pTi) • ••Fjv(ptjv) ! 



(Al) 



where N denotes the total number of outgoing particles. 
The Fj (px j ) are normalized to unity, and we further 
assume that they satisfy average momentum conservation 

(Pti H h Pt w ) = 0. 

Here, we also take into account exact transverse mo- 
mentum conservation, which is expressed by means of an 
additional Dirac constraint Ml : 



J 



dP 



dp Tl ■ --dp TN 



nN (p T ) Fi(pxi) • ••-Fjv(ptjv) 5 (Pti 



PTjv) 



(A2) 
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Note that Fj(pTj) no longer strictly represents the 
single-particle momentum distribution. The latter is now 
given by pi: 



dP 
dpT 7 



(A3) 



N(P 2 T ) )' 

We are going to compute the correlated probability 



distribution between the transverse momentum of an ar- 
bitrary particle (which we choose to be the one labeled 
N, Ptn) an d the flow vector Q. To compute this prob- 
ability distribution dP/dpTpfdQ, we int roduce the flow 
vector Q defined by Eq. (|J) in Eq. (A2) by means of a 
Dirac constraint, and integrate over the N — 1 remaining 
particle momenta: 



dP 



nN(p 2 T )F N {p TN ) /fll Fj(PTj)dp T )j 5{p Tl 



Ptat)<5 (wi 



+ w w - Q) , 



(A4) 



with Wj = WjUj = {wj/p Tj )p Tj . 

Using the Fourier representation of the Dirac distributions, this equation becomes 

JV-l 



dP 



dp T N dQ 



= nN(p* T )F N (p TN )J-±±^J (0^ 1Q IT (/ ' AP^-r^' dp T 



(A5) 



We have assumed that the particle with momentum Ptat is not involved in the determination of the reaction plane, 
i.e., that Wat = 0, as one usually does to avoid autocorrelations 0. 

Since N is large, the integrals over k and 1 can be evaluated by means of a saddle-point approximation, i.e., only 
small values of k and 1 contribute. Expanding to second order in k and 1 and re-exponentiating, one obtains 



J F(p T )e^ +l ^dp T ~ exp (ik • <p T ) + il • (w) - ^ (p%) - j (w 2 ) 



kl 



(wpj 



(A6) 



r 



where we have dropped the index j for simplicity. In 
deriving this equation, we have assumed that directed 
flow is small for most particles, \v%\ <C 1, so that terms of 
second order in the flow can be neglected, i.e. (pt) <C 
(Pt)) e t c - Note that the correlation between Ptn and 
Q comes from th e cr oss term, prop ortional to k ■ 1. 

Inserting Eq. ( |A6| ) into Eq. (A5), the remaining inte- 
grals over 1 and k are Gaussian, and can thus be evaluated 
analytically. After some algebra, and dropping the index 
N for sake of brevity, one finally obtains the result in the 
form 



dP 



dp T dQ 

l + 2c 



dP dP 

dpT dQ 
Pt (pt) 
Pt 



Q-(Q> 



(A7) 



In this equation, dP/ dp? denotes the single-particle mo- 
mentum distribution [obtained by integrating Eq. (A7) 
over Q], given by Eq. (A3), and dP/dQ the distribution 



of the flow vector Q [obtained by integrating Eq. (A7) 
over px] , given by Eq. (0) . The dimensionless correlation 
strength c is given by 



Pt 



f 



(A8) 



where / is given by Eq. (19). Finally, the width of the 
distribution of Q is 



N(w 2 )(l-f) 



(A9) 



In this expression, (1 — / 2 ) is the correction to the width 
due to momentum conservation. 

The above expressions of c and a have been obtained 
by a "brute force" calculation, but o ne c an easily recover 
them directly from the distribution (A7), and the corre- 
lation due to momentum conservation, Eq. (^), which we 
rewrite in the form 



(Ui • u 2 ) - (Ui) • (u 2 ) = 



PT1PT2 
'Nip 2 ,)' 



(A10) 



From the distribution of Q, Eq. (|7]), one obtains a 2 = 

(Q 2 ) - (Q) 2 - Using the definition of Q, E q (|) , and 
Eq. (AlC), it is a simple exercise to check Eq. (|A9|). Sim- 
ilarly, Eq. ( A7) allows to express the correlation strength 
c in terms of simple averages: c = ((pt • Q) — (Pt) ■ 
(Q)) / (ctPt) j which can also be easily calculated using 
Eqs. (|) and flAlOQ , 

As explained in Sec.|v|, there are two ways of elimi- 
nating the correlation in Eq. (A7). The first possibility 
is to shift the flow vector for each particle. Making the 
change of variables Q' = Q — cupx /pt in the distribution 
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Eqs. (0) and expanding to first order in the correlation 
c, one easily checks that the correlated term in Eq. (A7) 
cancels, so that the distributions of Q' and px are un- 
corre late d. Replacing c and a by their expressions (A9) 
and (A8), one recovers Eq. (|l5|). 

Instead of redefining the flow vector, one can compute 
the azimuthal correlation between the particle an d th e 
flow vector using the correlated distribution Eq. (A7). 
The first term in the right-hand side gives the correlation 
due to flow, while the second term corresponds to the 
correlation due to momentum conservation, that is, the 
term we wish to calculate. This term gives the following 
contribution: 



(cos(0 - $)) 



mom. cons. 



((Q)-(cosA$) e x -(Q)) 



(All) 



The resolution (cosA$) is given by Eq. (^), and (Q) 
is obtained by a similar calculation from the Gaussian 
distribution (Q): 

(0) = * ^ [(1 + X 2 ) Vx72) + X W/2)] . 

(A12) 



The terms proportional to 7i(x 2 /2) cancel, and one fi- 
nally obtains 



(cos(0 - $)) 



mom. cons. 



c ^ e -xV 2/o(x 2 /2) . (A13) 



In the limit \ <ti 1 (poor resolution), this equals approx- 
imately Cy/ir/2, while for \ ^> 1 (good resolution), the 
correction decreases like e/(2%). 



APPENDIX B: CORRELATION BETWEEN TWO 
SUBEVENTS 



We shall now compute the azimuthal correlation be- 
tween two subevents, which we denote Q a and Qb, in 
the presence of flow and momentum conservation. For 
simplicity, we assume throughout the following discus- 
sion that both subevents are equivalent, and moreover 
that they are of equal multiplicities. Since both contain 
a large number of particles, we may directly apply the 
central limit theorem. The normalized correlated distri- 
bution is thus Gaussian: 



dP 



C?Q a C?Qb 



M(i-c 2 ) 



exp 



(Q a - (Qa)) 2 + (Qb - (Qb)) 2 - 2C(Q a - (Q a )) ■ (Q b - (Q b )) 

<?l^-c 2 ) 



(Bl) 



In this equation, C represents the strength of the "non- 
flow" correlation (in particular, that due to momen- 
tum conservation) between the subevents. For C = 0, 



subevents are independent, and the distribution (Bl) fac- 
torizes into the product of two distributions of the type 
(9). In the more general case C^O, one recovers Eq. (|) 
by integrating Eq. (Bl) over Qb- 

When correlations due to momentum conservation are 
taken into account, the width cr 2 ub = (Q su b 2 ) — (Q, 



sub/ 

with the 



is given by a formula analogous to Eq. (A9) 
important difference that TV and f 2 are replaced by N/2 
and / 2 /2, since we use only half the particles to construct 
the subevents: 



sub 



Z! 

2 



(B2) 



Note that cr S ub is not simply equal to a/V2, because of 
the term arising from correlations. 

Let us now calculate the correlation strength C. It can 
be obtained by taking the following average value with 
the probability distribution (p|) M, |27j]: 



C 



((Qa - (Qa» ■ (Qb - (Qb»> 



(Qa ■ Qb) - (Qa) ■ (Qb) 



(B3) 



Using the definitions of Q a and Q b and the correlation 
from momentum conservation, Eq. (|A10 ), the numerator 
of Eq. (|B3|) is equal to -(N/2) 2 (wp T ) 2 / 'N {p 2 T ) , where 
we have neglected term s of second order in flow. With 
the width given by Eq. ([B2]), we finally obtain 



C 



(wp T y 



-/ 2 



2(p|; 



(w*) (1 - / 2 /2) 



r- 



(B4) 



where / is given by Eq. (|19|) . 

As explained in Sec. TlU the resolution parameter 
Xsub = (Qsub) • e x /cr su b of the subevent is obtained from 
the distribution of the relative angle between Q a and Qb- 
Unfortunately, there is no general analytic expression for 
this distribution. 

Useful expressions can be obtained in a few limiting 
cases. First, if the correlation strength C is much smaller 
than unity, the distribution (Bl) can be expanded to 
first order in C, which yields a distribution analogous 



to Eq. (A7). From this expression, one can in partic- 
ular evaluate the average cosine (cos($ a — $(,)), which 
leads to Eq. (pi]). Second, if the resolution is high 
(Xsub 2> 1), both Q a and Qb are close to their average 
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value (Q S ub)j so that their azimuthal angles with respect 
to the reaction plane are small, and one may approxi- 
mate $ a ~ Q a ■ e y /| (Q a ) |, and a similar expression for 
${,. In this limiting case, the distribution of <£> a and $t 
is a Gaussian, and the distribution of the relative angle 
is given by Eq. (p3|). 

Quite remarkably, one can derive an exact expression, 
valid for arbitrary Xsub and C , of the fraction of events 
for which the angle between subevents is larger than n/2, 
taking into account both flow and momentum conserva- 
tion. The condition |$ a — <£> b | > n/2 can equivalently be 
written Q a • Qb < 0, so that 



P(|$ -$ 6 | >tt/2) 



C?Q a rfQb 



dP 



rfQ a C?Qb 



0(-Qa-Qb). 

(B5) 



We then insert the Fourier representation of Heaviside's 
step function 6{x): 



9{-x) 



— lim 

2ir e^o+ 



+ 00 Jlikx 



-dk. 



IE 



(B6) 



With the probability distribution Eq. (Bl), the integrals 
on Q a and Qb are Gaussian and can be evaluated ana- 
lytically. One obtains 



P(|* n -*I,| >7T/2) = 



2tt 



+ OC 



dk 



(k + is) [1 - ik(l + C)] [1 + ik(l - C)] 
I 



exp 



2i X 2 k 



1 - ik(l + C) 



(B7) 



The only essential singularity of the integrand occurs in 
the lower half plane at k — — i/(X + C). Therefore the 
integral can be calculated by closing the contour on the 
upper half plane. The only pole is at k — — C), and 
the theorem of residues yields [p7fl 



P(|*«-$b| >?r/2) 



1-C 



"Xsub 



(B8) 



Inserting the correlation strength Eq. ( B4 ) in this expres- 
sion, one recovers Eq. (p2|). 
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